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Abstract. - We propose new polymer models for Monte Carlo simulation and apply them to 
a polymer chain confined in a relatively thin box which has both curved and flat sides, and 
show that either an ideal or an excluded- volume chain spends more time in the curved region 
than in the fiat region. The ratio of the probability of finding a chain in the curved region and 
in fiat region increases exponentially with increasing chain length. The results for ideal chains 
are quantitatively consistent with a previously published theory. We find that the same effect 
appears with excluded- volume chains and a similar scaling relation can be applied to them up 
to a certain length of the polymer. 



Introduction. - The structure Q, ||, j| and dynamics [||, [||, J?]] of polymer chains 
confined within a narrow space such as in a slit (Fig. |l|) are important for practical applications, 
e.g. filtration, gel permeation chromatography and oil recovery, and have attracted much 
interest. || The effect of confinement on polymer chains was calculated theoretically in well- 
defined geometries (a slit and a tube) and was simulated by a lattice model |J and an off-lattice 
model. jfij 

Recently, polymer confinement in curved boxes (a cylindrical shell, an undulated slit, etc) 
was studied, jnj], Such studies can relate to not only the behaviour of polymers confined 
in complex rigid geometries but also the behaviour of deformable bilayer membranes confining 
the polymer chains. 

Yaman et al . [fllf| showed that an ideal (Gaussian) polymer chain confined between cylindri- 
cal shells (Fig. |2|) has a lower free energy than one confined between two flat surfaces (Fig. p]), 
and predicted that polymer chains confined in bilayer membranes reduce the effective bending 
rigidity of the membranes and might induce spontaneous curvature in the system, e.g. leading 
to transitions from lamellar to bicontinuous phases. 
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Fig. 1. - flat box 




Fig. 2. - cylindrical shell 

Above mentioned studies of the confinement effects are limited to ideal chains. The ideal 
chain represents the polymers in solvents where the excluded- volume interactions are com- 
pensated for by monomer- monomer interactions. It can act as an unperturbed model for a 
perturbation calculation of the "real" chain in a good solvent characterized by excluded- volume 
interactions. However, the perturbation calculation is difficult for confined flexible polymers 
and few publications have addressed the excluded- volume effects. [|l3| 

In this paper, we investigate both curved surface effects and excluded- volume effects on a 
flexible polymer confined between two surfaces using Monte Carlo simulations and compare 
with an ideal chain confined in the same geometry. We suppose that the surfaces are purely 
repulsive (no trend towards adsorption), with and without excluded- volume interactions. In 
particular, we are interested in the polymer behaviour when the polymer chain has a large 
size relative to the gap between the confining surfaces. The effect is sensitive to the curvature 
of the surfaces, therefore, our off-lattice simulation has the strong advantage over a lattice 
simulation where the curved surface is described by a zig-zag boundary (steps of flat surface 
strips). 

First, we review the analytical study for ideal chains by Yaman et al. ]ll| Next, we propose 
new simulation models and check our algorithms in the free n-dimensional spaces to see if 
the models reproduce the known characteristics of flexible polymers. Then, we apply our 
simulation to ideal chains confined in a race-track box composed of curved and flat regions 
(Fig. ||) in order to compare with the analytical results by Yaman et al. Finally, we apply our 
simulation to excluded- volume chains confined to the race-track box. 



Analytical study for ideal chains by Yaman et al. - 



For a random walk of N steps with 
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Fig. 3. - race-track box 

endpoints at r and r', the propagator G(r, r', N) is defined as [|| 

{m ~ ? v G(r ' r '' JV) = 5(r " r ' )5{N) > (1) 

where Z is the length of one step. The ideal chain is modelled as a random walk in space, where I 
and N correspond to the characteristic length of a monomer and the degree of polymerizations, 
respectively. 

The partition function is a sum over all configurations: 

Z= V J drdr ' G ^ r ^ N )' ( 2 ) 
With the eigenfunction expansion: 

G(r, r', N) = Y, K(r)*n(r')e- NE ^ /6 e(N), (3) 

n 

where Q(N) is the step function, our problem of a polymer confined in a spherical annulus in 
arbitrary space dimension D with radius a and thickness d reduces to 

V 2 ^(r) = -EV(r), (4) 

under SO(D) symmetry and the boundary conditions: 

*(|r| =a) =*(|r| = a + d) = 0. (5) 

For the radial wave function: 

where e = Ed 2 , £ 2 = v 2 + (D — 2) 2 /2, and u = r/d, the boundary conditions are ^ u (v) = 
$>v(rj + 1) = 0, where rj = a/d. The general solution of Eq.@ is 

= AJ^y/hl) + BN^y/hl), (7) 

where J and A/" are the standard Bessel functions. The boundary condition determines the 
relation between e and rj. For spherical symmetry, only v = states contribute to the partition 
function. We label these s- wave states by n = 1, 2, . . ., and define the normalized radial wave 
function as R n (r) for v = 0. 
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For a spherical annulus with a large inner radius a and a small thickness d, E n can be 
obtained perturbationally: 



E n (a,d)=( K —) 



2 (D - !)(£> - 3) 
4d 2 



8n 2 7r 2 



(8) 



Note that D = 1 is a slit between two flat surfaces (Fig. |l|), D = 2 is a cylindrical shell (Fig. [|), 
D = 3 is a spherical shell. The curvature lowers the spherically symmetric energy spectrum 
for 1< D < 3. 

The partition function becomes 



1 CO 



-NE n (a,d)l 2 /6 



n=l 



y drr (jD - 1} ^ n (r) 



To indicate explicitly the dependence of the energies on a and d, we introduce B n (a,d): 

|2 



B n {a,d) = Sd _ i 



{{a + d) D -a D } 



(9) 



(10) 



where is the surface area of the unit sphere in D-dimensional space. The free energy is 
then given by 



F = -kT In 



^B n (a,d)e~ 



■E n (a,d)Nl 2 /6 



n=l 



(ii) 



In this paper, we are interested in the case of a cylindrical shell (D = 2) in three dimensions, 
for a polymer size R g ~ \ /f Nl larger than the shell thickness d. Then, the leading term of the 
free energy difference between the cylindrical shell and the flat box (D = 1) is 



Nl 2 I 2 
AF ~ kT(E 1 \ D=2 - E^) — = -kTN—. 



(12) 



From this result, Yaman et al. predicted that confining an ideal polymer chain into the 
bilayer membrane lowers the effective bending rigidity of the membrane and might induce 
spontaneous symmetry breaking of the membrane. Another interpretation of this result is 
that a polymer chain confined in the rigid shell is attracted into curved region, and that the 
ratio of probability finding a polymer chain in equivalent volume of curved and flat regions 
should be 



C 



-AF/kT 



exp 



/7V7 2 \ 
V24a 2 J ' 



(13) 



Here, we assume that the size of the confined polymer is smaller than the characteristic length 
of the curved and flat sections. 



Models and method. - Our polymer simulations are off-lattice beads and bonds models 
with N + 1 beads and TV bonds; they are relaxed versions of the random flight model with 
variable bond lengths. We call this the variable range bonding model. The models consist 
of the following conditions: (1) Connectivity: \ri — < 1.0, where = (r^i,r^2, • • • ,?m) 

stands for the position of the i-th bead. (2) Excluded- volume: \ri — Vj\ > 2r e , for all z, j. r e 
is the radius of the excluded volume of the beads. (3) Geometric constraint: the beads are in 
a confining geometry, here a race-track box (Fig. ph. Note that we impose the last condition 
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only on beads, but not on bonds. Thus, we even allow a bond that cuts through the curved 
wall as long as two beads attached to the bond satisfy the geometric constraint. We consider 
two models: Model-I which requires (1) and (3), and Model-E which requires (l)-(3). They 
correspond to the ideal polymer and the excluded- volume polymer respectively. 

Our Monte Carlo (MC) algorithm is the following, i) Select a bead at random, ii) Generate 
a trial jump: r — > r + A. Each component of the jump vector A M (/i = 1,2, • • • , n) is 
independently determined by the Gaussian distribution with zero mean value and a standard 
deviation of 0.3 for Model-I and 0.15 for Model-E. The Gaussian distribution is generated by 
the Box-Muller method. |Q iii) Check if the new position satisfies above mentioned conditions: 
(1) and (3) for Model-I, (l)-(3) for Model-E; if it does, accept the jump; if not, abandon the 
jump, iv) Return to i). For a simulation with TV + 1 beads, TV + 1 cycles comprise one MC 
step. 

The reason to halve the standard deviation of the Gaussian distribution of the jump distance 
for Model-E is to increase the acceptance ratio of the bead jump although the bead in Model-E 
is restricted by the excluded-volume interactions compared with Model-I. 

For Model-I in free n-dimensional space, the bond distribution is expected to be uniform 
inside the n-dimensional sphere with radius 1. Therefore, the bond length distribution is 
described by the distribution function P(r) of bond length r for < r < 1 as 

P(r) =nr n ~\ (14) 

and 

f p( r )^ r = i (15) 
Jo 

Thus, the mean bond lengths are analytically calculated as 

/ r 2 P(r)dr = n/(n + 2), (16) 
Jo 

and the chain should be a Gaussian chain with bond length ^n/(n + 2). The mean square 
end-to-end distance is analytically calculated as (R%) = 7Vn/(n + 2). For 3-dimensions, 
the mean square bond length becomes (/ 2 ) = 0.6. To evaluate the mean square end-to-end 
distances (R%) by simulations, we made ten runs of simulation with 4 x 10 6 MC steps for 
each length up to 200 beads. The simulation results are linearly fit by: (R%) = 0.3347V in 
ID, 0.4927V in 2D, 0.5947V in 3D, which implies that (/ 2 ) = 0.594 for 3-dimensions. They are 
consistent with the analytical calculations. 

For Model-E in the free n-dimensional space, we also carried out a series of simulations 
with 4 x 10 6 MC steps. We plot the results for free 3-dimensional space in Fig. 4. The fitting 
exponent v = 0.601 matches the best theoretical evaluation v = 0.588 with 2 % error. The 
square average of the bond length is obtained as (/ 2 ) = 0.6742. 

Our race-track geometry (Fig. ||) consists of two flat sections and two curved sections. The 
width of the track is d = 0.5, the radius of the inside wall of the curved sections is a = 4.75, 
and the length of the flat sections is L = 15.71. The volume of the curved region and the flat 
region are the same. 

For race-track simulations of Model-E, we use the same geometry as used for the Model-I, 
except that the width of the race-track d = 0.5 is defined relative to the center of the beads. 
When the width is defined relative to the surface of the beads, it is d + 2r e = 1.0 and the 
radius of the inside wall is a — r e = 4.5 where r e = 0.25 is the radius of the beads. 

To avoid localization on one side of the track, the polymer is relocated, every 10 6 MC steps, 
to the position of 180° rotation with respect to the center point of the track. The position data 
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Fig. 4. - Log-log plot of the end-to-end distance of Model-E in free 3D space. The fitting curve 
is i?oc TV - 601 . 
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Fig. 5. - Results of Model-I. The line is drawn according to the theoretical curve with param- 
eters a = 4.75, (Z 2 ) = 0.6: C = exp(7V/902.5). 



of beads are gathered every 1000 MC steps. Such data give us an idea of a density profile of 
the polymer in the restricted geometry. We divide the number of beads recorded in the curved 
region by the number of beads recorded in the flat region to obtain the probability ratio C. 

The advantages of our polymer models are 1) the movement of polymer in the narrow space 
is facilitated due to variable bond lengths and variable range bead jumping, and 2) the high 
probability of polymer approaching to the wall compared to models with fixed bond length 
and fixed range bead jumping is convenient to observe the effects in proximity to the curved 
wall. 

Results and discussions. - For each model and each polymer length, we carried out 5 runs 
of simulation with 4 x 10 8 MC steps after an initial randomized stage (10 6 MC steps). 

For Model-I, the race-track simulation is performed for each chain length (number of 
beads) : TV + 1 = 50, 100, 150, 200, 250. The results are plotted by in Fig. 5. The 
geometry gives parameters: d = 0.5 and a = 4.75. The mean square of the bond length 
is analytically calculated as (/ 2 ) = 0.6. Then, the theoretical curve with the parameters is 
given by C = exp(7V/902.5). 

For Model-E, the race-track simulation is performed for each chain length: N + 1 = 25, 50, 
75, 100, 125. The results are plotted in Fig. 6. We try to fit the data of the excluded-volume 
chains by the theoretical curve developed by Yaman et al. for the Gaussian chain. The geometry 
gives parameters: d = 1.0 and a = 4.5. The mean square of the bond length is evaluated by 
our free space simulations as (/ 2 ) = 0.6742. Then, the theoretical curve with the parameters 
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Fig. 6. - Results of Model-E. The line is drawn according to the theoretical curve for a Gaussian 
chain with parameters a = 4.5, (/ 2 ) = 0.6742: C = exp(7V/720.9). 

is given by C = exp(7V/947.5). 

In Model-I, our simulation results (Fig. 5) fit the analytical formula developed by Yaman 
et al. for all the chain lengths. In Model-E, our results (Fig. 6) are well described by the same 
formula as far as the chain lengths are short, however, they show a deviation when the chain 
gets longer. 

The excluded- volume effects result in a elongated equilibrium polymer shape compared with 
the ideal chain. This is confirmed by the snap shots of our simulations. They show that the 
length of the excluded- volume chain with N + 1 = 75 confined in the track is comparable to 
the section length L ~ 7ra. If the chain is too long to be accommodated in a single section, 
our interpretation of the free energy as a probability is not justified. Beyond this length, a 
single curved section cannot accommodate the whole chain and the chain cannot squeeze into 
the section due to the excluded-volume effects. Therefore, a deviation starting beyond this 
length is consistent with our analysis. The probability ratio C should reach a maximum in 
the vicinity of this length. This is consistent with our data. 

Conclusions. - 

We present new polymer models for Monte Carlo simulation which work well for our purpose. 
In particular, they have the advantage of efficient changes of the polymer conformation in 
restricted geometries due to the variable range character of bond lengths and of jump distances. 

We have investigated the curvature effect on a flexible polymer confined in a thin box. 
We have plotted the ratio of probability C finding a polymer between in a curved region 
and in a flat region as a function of chain length N. For the ideal chains, our results agree 
quantitatively with the prediction of the analytic theory. We also observe a similar curvature 
effect on excluded-volume chains whose lengths shorter than the length of the curved section 
of the box. 

Our simulations only cover the linear region of C as a function of chain length TV, due to the 
limitation of our computer power for the case of model-I and due to restriction of geometrical 
scale of confinement for the case of model-E. However, our results may still be useful under 
present circumstances where there is no experimental data for either ideal or excluded-volume 
chains and there is no analytical theory for excluded-volume chains. 

*** 

The authors thank F. Wittel, P. Pincus and K. Kremer for useful suggestions. Part of this 
work was done during the stay of one of the authors (Y.Y.S.) at MRL UCSB. 
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